Interaction induced staggered spin-orbit order in two-dimensional electron gas 
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We propose and formulate an interaction induced staggered spin-orbit order as a new emergent phase of 
two-dimensional Fermi gases. We show that when some form of inherent spin- splitting via Rashba-type spin- 
orbit coupling renders two helical Fermi surfaces to become significantly 'nested', a Fermi surface instability 
arises. To lift this degeneracy, a spontaneous symmetry breaking spin-orbit density wave develops, causing 
a surprisingly large quasiparticle gapping with chiral electronic states. Since the staggered spin-orbit order 
is associated with a condensation energy, quantified by the gap value, destroying such spin-orbit interaction 
costs sufficiently large perturbation field or temperature or de-phasing time. BiAg2 surface state is shown to be 
a representative system for realizing such novel spin-orbit interaction with tunable and large strength, and the 
spin- splitting is decoupled from charge excitations. These functional properties are relevant for spin-electronics, 
spin-caloritronics, and spin-Hall effect applications. 



PACS numbers: 71.70.Ej, 73.20.-r, 75.25.Dk, 79.60.Bm 

In semiconductor heterostructures, a charge particle mov- 
ing on a symmetry-breaking electric field experiences an 
effective 'anisotropic' magnetic field due to relativistic ef- 
fect, which couples to its spin. Such spin-orbit coupling 
(SOC), known as Rashba-dd or Dresselhaus-type[3 ] SOC, 
has proven to be a useful ingredient for realizing many 
physical concepts such as spin Hall effect l4l spin torque 
currentO |6), spin domain reversal,! 7] and exotic emer- 
gent superconducting, (8) and magnetic phases, [9 ] which are 
relevant for the applications of spin-electronics, 1 10] spin- 
caloritronics EH and quantum information processing. Along 
with long spin lifetime, spin relaxation, and immunity to 
perturbations or disorders, a highly desired recipe for these 
applications is decoupling the spin current from associated 
charge flow. This is a difficult task as the electric field gen- 
erated either by the broken inversion symmetry and/ or by the 
magnetic field itself naturally generates charge current per- 
pendicular to the spin current, thus greatly hinders the spin 
transport .1121 1131 Despite the discovery of giant SOC in a 
large class of condensed matter sy stems 121 17411191 practical 
realization and device implementation of these concepts have 
thus posed challenging. ifTQl fTTTl 

In this work, we introduce a theoretical proposal for gener- 
ating and manipulating staggered spin-orbit entangled helical 
state via electronic interaction. A key ingredient for realiz- 
ing such state is to have some form of intrinsic SOC, whose 
strength is not important, prior to the inclusion of interac- 
tion. For such systems, the ground state is defined by more 
exotic quantum numbers such as total angular momentum (J 
for j — j-type SOC) or helical quantum number (y = ± for 
Rashba- or Dresselhaus-type SOC), rather than typical spin- 
, orbital-, or momentum alone. Therefore, if the interac- 
tion strength is less than the SOC strength, a typical spin- or 
orbital-density wave alone is prohibited to form. On the other 
hand, we demonstrate here that novel emergent phases of mat- 
ter can arise in the spin-orbit channel even without necessarily 
breaking the time-reversal symmetry. 

We formulate the aforementioned postulates on the basis of 
a Rashba-type SOC ground state; however the idea is general 



and can be extended to other systems. Rashba SOC splits the 
non-interacting FS into two concentric helical Fermi pockets, 
and thereby a FS 'hot-spot' Q develops where degeneracy in- 
duces FS nesting between the opposite helical states. As a 
result of such instability, a translational symmetry breaking 
spontaneous collective ordering of helical degree of freedom 
develops, causing a giant SODW. The resulting SODW spa- 
tially modulates with a periodicity determined by the 'hot- 
spot' wavevector Q. The SODW renders gapping in the quasi- 
particle states, and the corresponding gap energy A thermo- 
dynamically represents its strength. Unlike in topological in- 
sulators or in other SOC systems where a spin-degenerate 
point exists with zero gapping, here the finite gap A protects 
the SODW phase from external perturbations and spin de- 
phasing. This means the spin- splitting survives up to a finite 
strength of Zeemen energy Ez (comparable to A), determin- 
ing the critical value of a time-reversal breaking perturbation 
such as magnetic field, B c , and spin de-phasing time r s . As 
a proposal, we show that a feedback effect of this order is 
the presence of a spin-orbit collective mode, which physically 
represents the interaction between two electrons with opposite 
spins and orbitals, or spin-orbit entangled quantum numbers, 
and that it can be detected via two-particle probes such as po- 
larized inelastic neutron scattering. 

Lattice model for Rashba coupling:-We start with a system 
of two-component Fermi gas in the presence of Rashba-type 
SOC. In the single-particle description of the system, the non- 
interacting FS is spin-polarized on the two-dimensional (2D) 
momentum space. For some systems, as in the surface state 
of BiAg2 deposited on various substrates, l2lfT9l the FSs can 
yield a shape to generate dominant FS nesting, and hence an 
unstable one-body ground state. If two momenta across the 
nesting 'hot-spot' is connected by spin flips, it cannot sup- 
port a charge density wave scenario. Furthermore, if the nest- 
ing commences between different segments of the same band, 
a spin-density wave may arise if the interaction strength can 
overcome the SOC strength. On the other hand, in such case 
a SODW can arise via inter-helical FS nesting instability even 
if the interaction strength is lower than the SOC strength. 
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FIG. 1. (Color online) (a) Non-interacting bands split by single- 
fermion Rashba-type SOC plotted along (lOO)-direction. The blue 
solid line is the coherent eigenstate, while the background is the 
associated single-particle spin-resolved spectral weight. The spin- 
polarization is depicted by a red (down- spin) to black (up- spin) gra- 
dient colormap. Along (100), the Rashba term vanishes, and thus the 
spin-resolved spectral weight gradient is absent here. The counter 
polarization of the spin texture along (110)- and (llO)-directions 
is demonstrated in (b) and (c), respectively, (d) Corresponding 
FS map is plotted in the same color scale, (e) Real part of non- 
interacting susceptibility at zero excitation energy reveals the devel- 
opment of paramount inter-helical state nesting at the 'hot-spot' vec- 
tor Q - 0.115(7r,7r). 



The non-interacting Hamiltonian in the two component 
fermion fields ^ k [^k^^k,]] T as H (k) ^Jj^fc 1 - 
^olr{ct x sin k y — a y sin k x )]ipk- Here is the free-fermion 
dispersion term, modeled by nearest-neighbor electronic hop- 
ping t as (fc = — 2t[cosk x + cos k y ] — E F , where E F 
is the chemical potential. The second term is the 2D lat- 
tice generalization of the standard Rashba SOC term H so = 
— zor(<x x k) z , with & being the Pauli matrices and (Xr being 
the Rashba SOC strength. The helical dispersion spectrum 
of H is E^ = £ fc ± a R [sin 2 k x + sin 2 k y ] 1/2 . The pa- 
rameters are obtained by fitting to the experimental dispersion 
from Ref.0 as listed in Ref.EUand the computed dispersion is 
plotted in Figs.[TJa),[TJb),[]J c ) a i on g different high-symmetry 
lines, and the corresponding FS is given in Fig.[TJd). 

FS instability: -We now investigate the FS instabilities of the 
system by evaluating the bare susceptibility \ m the particle- 
hole channel. \ is computed by convolving the single-particle 
Green's function G(k,iuj n ) = (iu n — ifo(fc)) -1 , yielding 
X(Q, Wm) = J2k,n iu n )G(k+q, iu n +ip m ), where iu n 



and ip m are the fermionic and bosonic Matsubara frequencies, 
respectively. The result plotted in the 2D q space at zero en- 
ergy (the real frequency is obtained by taking analytical con- 
tinuation from the Matsubara frequency) in Fig.[TJe) exposes 
that the nesting at Q = 0.115(7r,7r) is dominant. Q con- 
nects two momenta lying on different helical bands as shown 
in Fig.[TJd). Due to the definite chirality of the FS, different 
orientations of Q = 0.115(±7r, ±tt) vector are decoupled, 
and thus are exclusively included in the Hamiltonian. 

Spin-orbit density wav£.*-Based on these results, we now 
desire to write down a two-body interacting Hamiltonian 
in the reduced Brillouin zone in which the Nambu-Gor'kov 
spinor becomes ^ k [^k^^ki^k+Q^^k+Q,]} 7 • In 
this notation the interaction in the singlet-channel can be 
characterized by a contact interaction parameter g and cor- 
respondingly, ilj(fc) = g^t^k,&t + Q^k+Q,i> In or- 
der to reduce the two-body problem into ordered Fermionic 
ground state, we decouple the interaction term Hj by intro- 
ducing a auxiliary spin-orbit field A(fe) = gipk+Q v\ G z ® 
o-x]w'tpky [note that ip k ,v and ifrk+Qy belong to differ- 
ent helical states, having two Pauli matrices a z and a x for 
spin and momentum, and ® represents a tensor product be- 
tween them]. Employing mean-field approximation to the 
spin-orbit field A(fc), we obtain the total Hamiltonian as 
H = [H (k) + A(k)a z ® a x + h.c] which leads to 
the excitation spectrum as 

E£" = 5-+m[(^) 2 + A 2 ] 1 /2, 

with S£(fe) = ±££ )/2, (1) 

where v = — v = ± is the helical index due to SOC, and 
fi = ± is the split band index due to translational sym- 
metry breaking. We evaluate the order parameter A self- 
consistency as a function of temperature to a given contact 
potential. See supplementary material (SM) l22l for techni- 
cal details. The obtained temperature dependence and critical 
temperature demonstrate the spontaneous behavior of mean- 
field gap opening. 

Quasiparticle gapping: -Fig. [2] gives our main result in 
which the nature of the quasiparticle gap opening due to 
SODW is demonstrated and compared with angle-resolved 
photoemission spectroscopy (ARPES) data.|2] The blue lines 
are the coherent quasiparticle bands, plotted on top of the 
spin-resolved spectral weight [red (spin down) to black (spin 
up) colormap]. In the SODW state, two helical states split into 
several sub-bands, and a gap opens at the energy where two 
opposite helical states are connected by the Q vector. Along 
(lOO)-direction, the gap opening occurs in the filled state. The 
nature of gap opening and multiple number of shadow bands 
are in detailed agreement with recent ARPES data on the sur- 
face state of BiAg 2 alloy deposited on the monolayres (MLs) 
of Ag/Au(lll) heterostructure. l2l \W\ Similarly, along the di- 
agonal direction, A appears slightly above the Fermi level, 
and tiny hole-like pockets develop, as visible on the FS map in 
Fig.[2je). The constant energy surface map at an energy E=- 
110 meV illustrates how the main band and the folded band 
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FIG. 2. (Color online) Quasiparticle gapping due to SODW. (a) Electronic structure and associated spin-texture are plotted along (100)- 
direction. (b) Same as (a) but along (110) axis. The gap opening occurs below and above the Fermi level in both cases, while the spin- 
degeneracy at T-point is intact. Three horizontal arrows dictate three energy scales where the constant energy maps are presented in (d), 
(e), and (f). (c) The ARPES result of the gap opening in the BiAg2 surface state (obtained via deposition of the sample on varying number 
of monolayers of Ag/Au hetero structure), taken from Ref. Q. The present theory, replotted from (a) on top of the experimental data with 
red dashed line, is in detailed agreement with the ARPES data. The number of sub-bands developed in the SODW is clearly evident in the 
experimental data as well, demonstrating further the interaction origin of the gap. Results are presented in the unfolded Brillouin zone to 
facilitate comparison with experimental data. 



are nested along (100) -direction (in the particle-hole channel) 
where the quasiparticle gapping occurs, see Fig. [2jd). Sim- 
ilarly, at an energy E=62 meV above the Fermi level, the 
quasiparticle state is fully gapped along the diagonal direc- 
tion; Fig. gf). 

Tunable spin-orbit order:- Interestingly, the ARPES dataQ 
also reveals that the measured gap value A varies upon chang- 
ing the thickness of the Ag film, with a clear evidence of 
Kramers degenerate point at T-momentum for all cases. These 
observations indicate that the systems undergoes quantum 
phase transition with a spontaneous symmetry breaking in the 
spin-orbit channel, whereas the time-reversal symmetry is not 
explicitly broken. It is interesting to note that a similar broken 
symmetry state other than time-reversal symmetry is obtained 
on the surface state of topological insulators, l23l l24ll and also 
proposed to be responsible for the enigmatic 'hidden-order' 
state in heavy fermion URu2Si2,|25] despite having different 
forms of spin-orbit coupling quantum state in these systems. 

ARPES data plotted in Figs, gc), [Jf), [3ji) show that the 
gap decreases from A = 120 meV for ML2 to A = 80 meV 
for ML4 to A = 60 meV for ML 16. To theoretically explain 
this gap variation, we self-consistently tune the contact po- 
tential g, while keeping the Rashba term aR fixed. We find 
that the gap closing is associated with chemical potential vari- 
ation (beyond rigid band shift approximation), and the gapped 



region moves towards the Fermi level. We also find that the 
critical temperature below which the SODW sets in decreases 
with decreasing g\ see SM l22ll . Subsequently, the hole pock- 
ets along the diagonal directions grow in size. The large he- 
lical FSs with tunable area will be of considerable value for 
generating and detecting spin-resolved transports. In what fol- 
lows, the new interaction induced spin-orbit effect is tunable 
and large, even when the intrinsic SOC is fixed and weak. 

Spin-orbit correlation function and emergent collective 
mode: -Finally, to unravel the mechanism of SODW order, we 
deduce the emergence of associated collective excitations in 
the spin-orbit channel. The general form of the spin-orbit po- 
larization vector is defined as J(q,r) = Y^k^k v {t)[<j z ® 
v^w'^k+qyijY and correspondingly its bare correlator 
Xo(Qi t ) = (T T J(q,T)J(—q,0)), where r is the imaginary 
time and the operator T T denotes standard time-ordering be- 
tween the fermionic fields. The result is presented as a func- 
tion of excitation energy after performing Fourier transforma- 
tion, and including many-body correction within the treatment 
of random-phase approximation (RPA): x — Xo/0- ~ Uxo), 
where U is the interaction matrix defined in the SM[ 22 |. 

The corresponding result of the imaginary part of RPA sus- 
ceptibility x is given in Fig. 4. Along the diagonal direction, a 
robust collective mode is visible at Q ~ 0.115(7r, 7r) with its 
energy proportional to the gap value. The result is presented 
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FIG. 3. (Color online) Tunable spin-orbit order, (a) Quasiparti- 
cle dispersion along (lOO)-direction for self-consistently evaluated 
gap A=120 meV. (b) Corresponding FS map. (c) Comparison with 
ARPES data for the case of ML2.0 (d)-(f) Same as (a)-(c) but 
for a smaller gap of A = 80 meV and comparison with data for 
ML4 (reproduced from Fig. [2). The result for even a smaller gap of 
A = 60 meV is compared with experimental data for ML 16 in (g)- 
(i). As the gap gradually decreases with increasing number of mono- 
layers of Ag/Au heterostructure, the location of the gap approaches 
the Fermi level, suggesting that the interaction induced spin-order ef- 
fect is tunable, irrespective of a constant value of the non-interacting 
Rashba-coupling strength. 



for the case of ML2, however, the u res ~ A relation is re- 
produced for other cases (in the bare susceptibility level, the 
peak appears at the gap energy, however, the many body RPA 
correction shifts the mode energy to a slightly lower value). A 
downward dispersion of the mode is also visible centering Q, 
and extending to q — )> and q — » 0.124(7r,7r) with vanish- 
ing intensity. A second collective mode appears at an energy 
which is about half of the first mode energy due to multi-band 
effect. The first mode is a robust feature tied to the emergent 
spontaneous symmetry breaking SODW, whereas the inten- 
sity and energy of the second mode varies and its appearance 
is subject to the strength of interaction. The constant energy 
profiles at the first and second mode energies are depicted in 




[25 0.125-0.125 Q.125 
q=q y [n/a] (1/A) q x fr/a] (1/A) 

FIG. 4. (Color online) Collective spin-orbit mode, (a) Computed 
spectrum of imaginary part of RPA susceptibility in the spin-orbit 
channel in the SODW state plotted along (110) direction. The param- 
eters for this calculation is same as deduced for ML2 configuration 
before. A sharp peak in intensity at Q ~ 0.115(7r, tt) is clearly visi- 
ble around uj ~ A « 0.12 eV. A weak downward dispersion branch 
away from this mode can be marked. The second mode around A/2 
is a multiband feature, (b)-(c) Constant energy profiles of the same 
susceptibility at first and second mode energies, respectively. 



Figs. 4(b), (c). The bosonic spin-orbit excitation represents 
electron-electron interaction with simultaneous spin and or- 
bital flips, or an entangled helical index flip. [22] A polar- 
ized inelastic neutron scattering measurement, which directly 
probes the imaginary part of the susceptibility, will be able to 
detect this mode. We note that no such collective mode de- 
velops in the charge susceptibility, and thus confirming that 
the present spin-orbit interaction does not activate any charge 
excitation. 

Outlook and discussion: -The proposed spin-orbit order is a 
novel quantum phase of matter which can also arise in va- 
riety of other materials in which the non-interacting wave- 
function is defined by exotic quantum number such as heli- 
cal index, pseudospin, total angular momentum (typically in 
heavy-fermions and actinides) owing to SOC of various na- 
tures. Depending on the nature of the interaction and broken 
symmetry, the spin-orbit order can as well emerge as short 
range order in which other interacting phases such as quantum 
spin-Hall effect, spin-orbit nematic phase may exist. Here, we 
take a 2D electron gas as a representative example of how a 
SODW can arise due to the FS instability in a Rashba-type 
SOC background. In this context, it is interesting to point out 
the recent experimental findings of quasiparticle gapping in 
the surface state of topological insulator due to quantum phase 
transition. Such gap opening in the absence of time-reversal 
symmetry breaking and without the destruction of bulk topo- 
logical properties, 1 23, 24 ] violates the conventional topologi- 
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cal paradigm and criterion. (26) We envisage that it is sugges- 
tive of an emergent time-reversal invariant spin-orbit order, 
whose detail is required to be formulated in future study. 

Some important advantages of the present spin-orbit or- 
der than the typical single-particle SOC can be noted. (1) 
In topological insulators or in Rashba systems, the functional 
use of spin-orbit effect is subject to counter propagating he- 
lical spin state which is topologically protected by Kramers 
spin-degeneracy with zero gap at the r-point.|26] Therefore 
the presence of any type of time-reversal breaking impurity 
or defects, with strength as small as infinitesimal value will 
destroy the protection, and thus barring the exciting useful- 
ness of topologically protected transport properties. On the 
other hand, the present spin-orbit order is not only protected 
by symmetry, but most importantly it is thermodynamically 
shielded with a condensation energy equal to the tunable and 
large gap. Therefore, it would require a sufficiently large 
value of magnetic field such that its associated Zeeman en- 
ergy Ez ~ 9HbB (g is 'g' -factor, hb is Bohr magneton)] can 
overcome the condensation gap energy. (2) Spin de-phasing 
time r s which is an important ingredient for the spintronics 
and quantum computing applications is considerably larger 
here (determined by Ez energy), and is tunable. (3) Another 
crucial benefit of the interaction induced spin-orbit effect is 
that, unlike in typical electric field or magnetic field induced 
SOC, it does not necessarily activate a charge flow, and thus 
will be highly valuable for solely generating spin-transport. 

Another experimental verification of the broken symme- 
try spin-orbit order is detection of spin Nernst effect. Since 
Nernst effect is sensitive to reconstructed FS, a manifestation 
of broken symmetry phase, El [28) spin-orbit order will gen- 
erate a spin-resolved thermal current which are detectable in 
recent days laboratory facilities. (29) 

The author is grateful to M. J. Graf, A. V. Balatsky, J.-X. 
Zhu, S. Raghu, P. Wolfle for many fruitful discussions. The 
work is supported by the U.S. DOE through the Office of Sci- 
ence (BES) and the LDRD Program and faciliated by NERSC 
computing allocation. 
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SUPPLEMENTARY MATERIAL 

In the supplementary mayerial we provide the detailed 
derivation of the spin-orbit order parameter. Due to the 
definite chirality of the FS, different orientations of Qi = 
0.115(±7r, ±7r) are decoupled and are exclusively included 
in the Hamiltonian. In this spirit we define the Nambu- 
Gor'kov spinor as = [^ fc>t , ijj k , h ^k+Q^^k+Qi^ ...] T , 
where i = 1 — 4. Since macroscopically the gap value is 
same for all four Qi vectors, the interaction potential is con- 
stant which gives the two-body spin-orbit term as Hj(k) = 

^Z)i^fc,t^fc,^fc+Q 4 ,t^+Qi^- In the functional P ath in- 
tegral formalism, the partition function of the Fermi gas is 

Z = J V^V^ exp (-S[tf t, *]) (fi = k B = 1 throughout), 

where the action is *] = $ dr E'/J^ t<9 r^ + #o + 

Hi], with (3 = 1/T and the prime over the summation repre- 
sents that the summation is performed in the reduced Brillouin 
zone. 

In order to reduce the above-derived action into an or- 
dered Fermionic problem, we decouple the four-field in- 
teraction term by introducing the auxiliary spin-orbit cou- 
pling field A(fc) = gip{+ Q ^[o- z <g> (J x ]vv>^ky. Em- 
ploying the Hubbard-Stratonovich transformation! 1] and 
integrating out the fermionic fields, we obtain Z — 
|PA t PAexp(-6' e// [A t , A]), where the effective spin- 
orbit order reads S e// [At,A] = — |Trln [G _1 (fc,r)] - 
Jo Efc dr\A(k)\ 2 /g. The inverse single-particle Green's 
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FIG. 5. Self-consistent gap values as a function of temperature for 
A(0)/|£| ~0.1 with g/\t\=l. The critical temperature for the spin- 
orbit density wave, T so , is estimated for this case to be around 800 K. 
The line depicts a fit of the function A (T)/ A (0) = (l-T/T so ) 35 , 
which is close to the BCS critical exponent of 0.5. Inset: The critical 
temperature T so (blue squares) and gap amplitude A(0) (red circles) 
varies linearly with the interaction potential strength g. The data are 
presented with respect to \t\. 



function in the fermionic Matsubara frequency is obtained as 



G~ 



iuj n l - H (k) 

a T A^ iujr> 1 



CTj-A 

■ffo(fe + Qi) 



(2) 



The other terms belong to the three values of Qi. At the mean- 
field level, A(fe) represents the gap parameter which is same 
for all values of Qi. In this case, the effective action fur- 
ther simplifies to S e ff[A^A] = — \ In [detG -1 (fc, iuo n )) — 
Nf3\A(k)\ 2 /g, where N is the size of the system. From 
det[G -1 (fc, E)} = 0, we derive the excitation spectrum as 
given in Eq. 1 in the main text. Finally, by using the relation 
for the thermodynamic potential Q = — 1/(3 hi Z, we obtain 
the self-consistent condition for the number operator and the 



gap function as: 

n t/i = Yl / ^-ImGn/22(fe,^^ + ^)/M, ( 3 ) 

k v J-oo *K 

A= 9l2 ^-vImG 12 (k,v,u + i5)f{Lu), (4) 

where S is infinitesimally small broadening parameter as f(u) 
is the Fermi-Dirac distribution function for the excitation en- 
ergy of uj. We solve Eqs.[3]and|4] self-consistently to obtain 
the value of quasiparticle gap to match the experiment. 

The values of all parameters are obtained by fitting the 
quasiparticle-spectrum and gap with the experimental value 
from Ref. 0. We find t = -1.216 eV, and a R = 1.25/\t\. 
The interaction potential is self-consistently evaluated to be 
g/\t |=1, 0.84, and 0.73 which gives the experimental gap val- 
ues of A =0.12, 0.08 and 0.06 eV for ML2, ML4 and ML 16, 
respectively. We also deduce the temperature dependence of 
the gap parameter as given in the supplementary Fig. [5] The 
result reproduces a mean-field like behavior with the critical 
exponent 0.35, which can be compared with the typical BCS 
value of 0.5. The critical temperature below which a spin- 
orbit order sets in relies on the contact potential (see inset to 
Fig.0. 

SPIN-ORBIT EXCITATION 

Finally, we expand on the deduction of the spin-orbit sus- 
ceptibility. The spin-orbit operator is defined in the main text 
as 

^k+qyij). (5) 

k 

Substituting the expression for spin-orbit polarization in the 
spin-orbit susceptibility formula, we obtain 



Xo(Q,r) = -l(r T J(q,T)J(-g,0)> = ~ £(T t V>£ i „(t)V*+, xWl>(0)Vv- g x(0)> 



k,k' 



(6) 



I 

Rewriting the above susceptibility in the reduced Brillouin zone (denoted by prime over the summation), we get 

I 



1 ' \ 

k,k' L i 

(0)})] 



(7) 
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(b) 



k, co, v 




k+q, co+n, v' 



FIG. 6. (a) Real-space view of spin-orbit density wave having two 
opposite spin-orbit states sitting in different sublattices. Arrows 
stand for spin while circles with different colors give different or- 
bitals. We chose a commensurate modulation here along the diagonal 
for illustration convenience, (b) Spin-orbit bubble which is responsi- 
ble for electron-electron interactions having opposite spin and orbital 
quantum numbers. The interaction vertex Y vv i — [a z <j x ]w>- 

Expanding the above four-fermionic fields into two ordered 
fields in the spin-orbit channel and employing the Wick's 
theorem, we can write down the above equation in terms 
of Green's functions as 

Xo{q,r) = - — ^^G ij (fc,zy,r)G ij (fc + g,z/ / ,-r), 



k ij 



where Gij are the components of the Green's function given 
in Eq. 1 above. Here we employ the momentum conserva- 
tion relation k' = k + q or k' = k + q + Qi. Each term 
in xo represent a polarization bubble as represented in supple- 
mentary Fig.[6jb) which is responsible for interaction between 
two electrons with opposite orbitals and spins, or with oppo- 
site helical index v. The explicit form of the susceptibility ten- 
sorisxf (q,i Pm ) = -1/N E fc , w ' <t>$(k)4i(k)<$(k+ 
q)^, (fc + g)(/(25£) - f(E^' +q )/(ico n -E»- E<£ +q ). Here 
and are the eigenvector and eigenvalue of the Hamil- 
tonian deduced in Eq. [2] with /i, // being the corresponding 
band indices, and i,j,k,l their orbital counterparts. Within 
RPA framework, the perturbation interaction term for the sus- 
ceptibility is H RPA = UYj V n v n v + U'J2v n v n v> where 
v = — v = ± is the helical index and U, and U' are the in- 
teraction strength between same and opposite helical states. 
Since the interaction here is Coulomb type which acts on 
charge particles, therefore, without any loss of generality, we 
take U = U' ~ g = 0.35 eV for ML2 configuration. It is in- 
teresting to note here that despite the same value of interaction 
between both helical states, the long-range order is prohib- 
ited along the spin-channel due to strong spin-orbit coupling 
(U « olr), but from symmetry argument, an order parame- 
ter in the spin-orbit channel ([/'-term) is allowed. 



(8) 
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